Lifetimes of image-potential states on copper surfaces 
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The lifetime of image states, which represent a key quan- 
tity to probe the coupling of surface electronic states with the 
solid substrate, have been recently determined for quantum 
numbers n < 6 on Cu(100) by using time-resolved two-photon 
photoemission in combination with the coherent excitation 
of several states (U. Hofer et al, Science 277, 1480 (1997)). 
We here report theoretical investigations of the lifetime of 
image states on copper surfaces. We evaluate the lifetimes 
from the knowledge of the self-energy of the excited quasi- 
particle, which we compute within the GW approximation 
of many-body theory. Single-particle wave functions are ob- 
tained by solving the Schrodinger equation with a realistic 
one- dimensional model potential, and the screened interac- 
tion is evaluated in the random-phase approximation (RPA). 
Our results are in good agreement with the experimentally 
determined decay times. 



In a variety of metal surfaces, the bulk band structure 
projected onto the surface presents a band gap near the 
vacuum level. Thus, given an electron located outside 
such a metal surface, it may be trapped in the vacuum 
well produced by the self-interaction of the electron with 
the polarization charge it induces in the surface . Far 
from the surface, into the vacuum, this potential well ap- 
proaches the long-range classical image potential, and the 
resulting quantized electronic states form a Rydberg-like 
series which converges towards the vacuum energy and is, 
in principle, resolvable ]3J. These so-called image states 
were first identified experimentally by inverse pho- 
toemission |5|J|, and the first high-resolution measure- 
ments of image states were performed by the use of two- 
photon photoemission (2PPE) [0-0. Combined with 
the use of ultrafast lasers, 2PPE has provided a powerful 
technique to directly probe, on a femtosecond time scale, 
electron dynamics of excited electrons in metals [ p]Jl2| 
and, in particular, image states |12-n|. 



Besides their well-defined physical properties, image- 
potential induced states are of general interest in sur- 
face science, because of their applicability to other areas 
of condensed-matter physics, providing, in particular, a 
very simple model to investigate the coupling of surface 
electronic states with the solid substrate directly in the 
time domain |lq] . The knowledge of this coupling is cru- 



cial to understand many electronically induced adsorbate 
reactions at metal surfaces fL7| |) and it is this coupling 
with the underlying substrate which governs the lifetime 
of the image states. 

By using time-resolved 2PPE the lifetime of image- 
potential states on copper surfaces has been recently de- 
termined [^2|dl5| . Lifetimes of electrons in the n = 1 im- 
age state on Cu(lll) have been investigated [^2|dl4[, and 
Hofer et al |l5| have used time-resolved 2PPE in combi- 
nation with the coherent excitation of several quantum 
states to investigate the ultrafast electron dynamics of 
image states on the (100) surface of copper, providing 
an accurate experimental determination of the lifetime 
of the first six image states on that surface. 

A quantitative evaluation of the lifetime of image states 
was first reported in Ref. [Tq| , within a many-body free- 
electron description of the metal surface and neglecting, 
therefore, any effects of the ion cores. Hydrogenic-like 
states with no penetration into the solid were used to 
describe the image-state wave functions, a step model 
potential was introduced to calculate the bulk final-state 
wave functions, simplified jellium models were used to 
approximate the screened Coulomb interaction, and the 
image-state lifetimes were evaluated, within the GW ap- 
proximation of many-body theory |l9|], from the knowl- 
edge of the electron self-energy. More realistic image- 
state wave functions were introduced in subsequent cal- 
culations |^o| , allowing for penetration into the crystal. 

In this letter we present the results of a calculation of 
the lifetimes of the first three image states on a Cu(100) 
surface and the first image state on Cu(lll), by go- 
ing beyond a free-electron description of the metal sur- 
face, and including, therefore, band structure effects. 
First, we evaluate the electronic wave functions by solv- 
ing the time-independent Schrodinger equation with a 
realistic one-dimensional model potential, we then use 
these wave functions to evaluate the screened Coulomb 
interaction within a well-defined many-body framework, 
the random-phase approximation (RPA) fl2lfl , and we fi- 
nally evaluate the lifetimes from the knowledge of the 
imaginary part of the electron self-energy of the excited 
quasiparticle, which we compute within the so-called GW 
approximation. 

As the image-state wave functions lie mainly in the vac- 
uum side of the metal surface, and the electron moves, 
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therefore, in a region with little potential variation paral- 
lel to the surface p2^ ] , we assume translational invariance 
in the plane of the surface, which is taken to be normal to 
the z axis, and evaluate the damping rate of an electron 
in the state (/>o(z)e ik r r ii with energy E = (E ) z + kjj/2 

(we use atomic units throughout, i.e., e 2 = Ti — m e = 1), 
from the knowledge of the two-dimensional Fourier trans- 
form of the electron self-energy S(z, z'\ kii , Eq), as follows 



r = — 



2 J dz J dzVo(z)ImS(z,z';k||,£;o)0o(z')- (!) 



In the GW approximation, only the first term of the 
expansion of the self-energy in the screened interaction 
is considered, and after replacing the Green function (G) 
by the zero order approximation (G°), one finds: 



ImE(z, z'; k|| , Eq 
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where the sum is extended over a complete set of final 
states </> / (z)e i(k ii +q ii ) ' r ii with energy E f = (E f ) z + (k|| + 
q||) 2 /2, E F is the Fermi energy, and W md (z, z'\ qy, E) is 
the induced interaction: 



^ lnd (z,z';q||,£) = J dzt J dz 2 V(z,z 1 ;t lll ) 

xx(zi,z 2 ;q\\,E)V(z 2 ,z / ;q\\). (3) 

Here, V{z,z'\q\\) and x( z > z '\ q\\, E) represent two- 
dimensional Fourier transforms of the bare Coulomb po- 
tential and the density-response function, respectively. 

In the RPA, x( z , z '\ E) satisfies the integral equa- 
tion 



x{z,z';q ]l ,E) = x {z,z';qi\,E) + J dz x J dz 2 

x x °(z, z\\ q|| , E)V(z u z 2 ; q\\)x( z 2, z '\ q\\ , E), (4) 

where x°( z i: z 2] q\\ , E) represents the density-response 
function for non-interacting electrons. An explicit ex- 
pression for x°( z ii z 2\ q\\ , E) can be found in Ref. [p3| , 
in terms of the eigenfunctions 4>i(z) of the one-electron 
effective hamiltonian. 

We compute the image-state wave function, </>o(z), and 
both the final states <j>f(z) and all the one-electron eigen- 
functions involved in the evaluation of the polarizabil- 
ity X ( z it z 2] q\\ )E), by solving the Schrodinger equation 
with the one-dimensional model potential recently sug- 
gested in Ref. [M, which approaches, far outside the 
surface, the classical image potential and describes, in- 
side the crystal, the self-consistent effective potential of 
density-functional theory. This one-dimensional poten- 
tial reproduces the width and position of the energy gap 
at the f point (k|| = 0) and, also, the binding energies of 
both the n = crystal induced surface state at f and the 



first (n = 1) image-potential induced state. This model 
potential reproduces first-principles calculations of the 
wave functions and binding energies of image states in 
Li(100) and also the average probability-density of 
the n = 1 image state on Cu(100) derived by Hulbert et 
al Q from a first-principles calculation (see Fig. 1). We 
note that the members of the Rydberg series on Cu(100) 
with quantum number n < 3 have binding energies of 
0.57, 0.18 and 0.08eV. The n = 1 probability-density 
has a maximum at 3.8A outside the crystal edge (z = 0), 
which we choose to be located half a lattice spacing be- 
yond the last atomic layer. For the fraction of the first 
image state overlap probability with the bulk [z < 0) 
crystal, we find p = 0.05. On the Cu(100) surface the 
first image state is closer to the center of the band gap 
than on the Cu(lll) surface. As a result, on the latter 
the penetration of the first image state into the crystal is 
larger, p = 0.22, and the probability-density has a maxi- 
mum closer from the surface, at 2. 3 A outside the crystal 
edge. The binding energy of this image state is 0.83eV. 

First of all, we focus on the evaluation of the damping 
rate of the n — 1 image state of Cu(100), and we set 
the wave vector of the image electron parallel to the sur- 
face, k|| , equal to zero. Coupling of the image state with 
the crystal occurs through the penetration of the image- 
state wave function into the solid and, also, through the 
evanescent tails of bulk states outside of the crystal. Ac- 
cordingly, we have calculated separately the various con- 
tributions to the damping rate by confining the integral 
in Eq. (1) to either bulk [z < 0) or vacuum (z > 0) 
coordinates, and we have obtained the results presented 
in Table I. The decaying rates in the bulk are expected 
to be larger than in the vacuum. Also, for final-state 
wave vectors parallel to the surface, qy , that are smaller 
than y/2(Eo ~ E g ) (E g represents the lower edge of the 
projected energy gap at the f point) there is a reduc- 
tion in the phase-space because of the presence of the 
projected band gap, and this results in a reduced decay 
rate in the vacuum for which vertical transitions (qy w 0) 
are expected to dominate. Thus, the coupling of image 
states with the crystal states occurring through the bulk 
penetration plays an important role in the determina- 
tion of the damping rate, though the penetration of the 
image-state wave function is small. The contribution to 
the damping rate coming from the interference between 
0^0 and z' >0 coordinates is comparable in magnitude 
and opposite in sign to both vacuum and bulk contribu- 
tions, this being a consequence of the behaviour of the 
imaginary part of the two-dimensional Fourier transform 
of the self-energy [ p5[ . The experimentally determined 
damping rate reported by Hofer et al |15| for the first 
image state on Cu(100), also presented in Table I, shows 
a reasonable agreement with our theoretical prediction, 
though our predicted linewidth is a little larger. 

Also exhibited in Fig. 1 is the square of the 
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parametrized hydrogenic-like wave function <f>(z) = 
(4a 3 ) 1 / 2 ze~ az , with the a parameter chosen so as to 
fit the peak position of our calculated first image-state 
probability-density of Cu(100). This hydrogenic-like 
wave function is found to be less localized around the 
maximum outside the surface, and replacement of our 
calculated image state wave function by this approxima- 
tion leads, therefore, to a damping rate (see Table I) that 
is too large. First, the use of this approximated initial 
wave function gives rise to a spurious contribution from 
a region just outside the surface, in which the decaying 
rates are noticeably larger than around the maximum, 
and we note that the damping rate is highly sensitive to 
the actual shape of the image-state wave function. Sec- 
ondly, both the bulk contribution and the contribution 
from the interference between z <0 and z' >0 coordinates, 
which we have found to be of crucial importance, are 
completely neglected within this model. 

In order to investigate the impact of band structure 
effects on the damping rate of image states, we have also 
performed a calculation in which the realistic final states 
4>f(z) considered above are replaced by the self-consistent 
jellium LDA one-electron wave functions of Lang and 
Kohn p^] , with and without the restriction that only 
final states with energy (Ef) z lying below the projected 
band gap are allowed. Both the image-state wave func- 
tion and all the one-electron eigenfunctions involved in 
the evaluation of the screened interaction are still ob- 
tained by introducing the one-dimensional potential of 
Ref. pij into the effective hamiltonian. The various con- 
tributions to the damping rate calculated in this way 
are presented in Table I. For > ^2(Eq — E g ) all fi- 
nal states with energy Ef < Eq lie below the projected 
band gap; thus, the bulk and interference contributions 
to the damping rate remain almost unaffected by this 
restriction. However, as the coupling of the image state 
with the crystal occurring through the tails of bulk states 
outside of the crystal is expected to be dominated by ver- 
tical transitions (qn ks 0), the vacuum contribution to the 
damping rate becomes dramatically smaller as final states 
lying within the projected band gap are not allowed. 

As far as the screened interaction is concerned, sys- 
tematic investigations of the role that this quantity plays 
in the coupling of image states with the solid have been 
performed |27]]. We find that simplified jellium models 
for the electronic response lead to unrealistic results for 
the lifetime, though the impact of the band structure in 
the evaluation of the screened interaction is not large. 

In a heuristic view, the bulk contribution to the de- 
cay rate for an image state might be approximated by 
the value of the decay rate for a bulk state at the same 
energy, times the fraction of the image state overlap prob- 
ability with the bulk crystal. We have calculated, within 
the GW-RPA, the decay rate for bulk states in a homo- 
geneous electron gas, thereby neglecting band structure 
effects in both the initial <f>o(z) and final states <fif( z ) and, 



also, in the polarizability. The damping rate for a bulk 
state at the energy of the first image state on Cu(100), 
times the penetration of this image state, results in a 
linewidth of 14meV pfi|] , well below the bulk contribu- 
tion to our full GW-RPA surface calculation presented 
in Table I. 

Finally, we consider the lifetime, r, of image-potential 
states on Cu(100) with quantum number n < 3, and 
the first image state on Cu(lll). As before, image- 
state wave functions, final states and all the eigenfunc- 
tions involved in the evaluation of the polarizability are 
obtained by solving the Schrodinger equation with the 
one-dimensional model potential of Ref. 0] , and we set 
k|| = 0. The results of these calculations are presented 
in Table II, together with the experimentally determined 
decay rates reported in Refs. (l^Jl^l and for Cu(lll) 
and Cu(100), respectively. We observe that both our 
calculated and the experimentally determined lifetimes 
of the image states on Cu(100) increase rapidly with n 
as ~ 1/rt 3 , in agreement with previous theoretical pre- 
dictions |jj ; this is a result of the weaker spatial overlap 
with the bulk that higher order image states present. As 
for the lifetime of the first image state on Cu(lll) we 
find that it is comparable to that of Cu(100), though 
this image state on Cu(lll) is located close to the top 
of the projected band gap and the wave function overlap 
with the bulk is, therefore, much larger. Also, the surface 
electronic structure of Cu(lll) supports the well-known 
ri = crystal-induced surface state, and we find that the 
decaying rate of the n = 1 image state to this intrin- 
sic surface state results in a linewidth of 16meV, which 
represents a 40% of the total linewidth (r _1 = 38meV). 
However, both the large bulk-state overlap and the ex- 
istence of the n = surface state are counterbalanced 
by the band gap extending on the Cu(lll) surface below 
the Fermi level and the available phase space becoming, 
therefore, highly restricted. 

In conclusion, we have reported theoretical calculations 
of the lifetime of image-potential states on copper sur- 
faces, by going beyond a free-electron description of the 
metal surface. We have analyzed the origin and mag- 
nitude of the various contributions to the quasiparticle 
damping. Our results indicate that contributions to the 
damping rate coming from the bulk and from the inter- 
ference between bulk and vacuum coordinates both play 
an important role in the determination of the lifetimes. 
We have also demonstrated that the vacuum contribu- 
tion to the damping rate is highly sensitive to the lo- 
cation of the projected band gap. Our results are in 
reasonably good agreement with the experimentally de- 
termined decay times reported by Hofer et al |l5| in the 
case of Cu(100), and our calculated lifetime of the first 
image state on Cu(lll) is in excellent agreement with 
the experimental results reported in Refs. jl3| and |3 . 
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just above the Fermi level (Ef) the damping rate is found 
to be (see J. J. Quinn and R. A. Ferrell, Phys. Rev. 
112, 812 (1958)) r" 1 = 0.00251rf /2 (£ - EpfeV' 1 , in 
linewidth units. This approximated formula gives for the 
damping rate for a bulk state at the energy of the first im- 
age state on Cu(100), times the penetration, a linewidth 
of 24meV. 

TABLE I. Damping rates, in linewidth units (meV), of the 
n — 1 image state on Cu(100). (a) Our full calculation, as 
described in the text, together with the experimental damp- 
ing rate of Ref. Jlq |. (b) The result of our calculation when 
our realistic image-state wave function is replaced by the hy- 
drogenic-like wave function described in the text, (c) and (d) 
The result of our calculation when our realistic final states are 
replaced by the self-consistent jellium LDA wave functions of 
Lang and Kohn JifJ, (c) with and (d) without the restric- 
tion that only final states with energy {Ef) z lying below the 
projected band gap are allowed. 





Bulk 


Vacuum 


Interference 


Total 


Experiment 


(a) 


24 


14 


-16 


22 


16.5 


(b) 




71 




71 




(c) 


19.5 


13 


-11 


21.5 




(d) 


21 


58 


-12 


67 





TABLE II. Calculated lifetimes, in femtoseconds, of the 
n < 3 image states on Cu(100) and the first image state on 
Cu(lll), together with the experimentally determined life- 
times of Refs. and (TH) for Cu(lll) and Cu(100), re- 
spectively. The damping rate of Eq. (1) is related to the 
lifetime, r, by r _1 r = la.u. = 660meVfs. 



Cu(100) 
Theory Experiment |l5 


Cu(lll) 
Theory Experiment fillJll 


n = 1 30 40 ± 6 


17.5 18 ±5 


li 




15 ±5 




n = 2 132 110 ±10 






n = 3 367 300 ± 15 







FIG. 1. Probability-density of the n = 1 image state on 
Cu(100) (solid line) and the hydrogenic-like approximation 
described in the text (dashed line) . Stars represent the prob- 
ability-density of the n = 1 image state, averaged parallel to 
the surface, as reported in Ref. [E2j. 
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